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Abstract. We investigate the nature of the critical behaviour of the random- 
anisotropy Heisenberg model (RAM), which describes a magnetic system with random 
uniaxial single-site anisotropy, such as some amorphous alloys of rare earths and 
transition metals. In particular, we consider the strong-anisotropy limit (SRAM), 
in which the Hamiltonian can be rewritten as the one of an Ising spin-glass model with 
correlated bond disorder: Ti = —JJ2{xy)-^xyO'x<^y, where jxy = Ux ■ Uy and is a 
random three-component unit vector. We performed Monte Carlo simulations of the 
SRAM on simple cubic lattices, up to i = 30, measuring correlation functions 
of the replica-replica overlap, which is the order parameter at a glass transition. 
The corresponding results show critical behaviour and finite-size scaling. They 
provide evidence of a finite-temperature continuous transition with critical exponents 
ijo = —0.24(4) and i^o = 2.4(6). These results arc close to the corresponding estimates 
that have been obtained in the usual Ising spin-glass model with uncorrelated bond 
disorder, suggesting that the two models belong to the same universality class. This 
is consistent with arguments that suggest that the disorder correlations present in the 
SRAM are irrelevant. 



PACS numbers: 75.50.Lk, 05.70.Jk, 75.40.Mg, 77.80.Bh 
1. Introduction 



Extensive theoretical and experimental work has been devoted to the study of 
amorphous alloys of rare earths and transition metals, for instance TbFe2 and YFe2. 
They are modeled by a Heisenberg Hamiltonian with random uniaxial single-site 
anisotropy defined on a simple cubic lattice, or, in short, by the random-anisotropy 
model (RAM) 

H = -J^Sx ■ Sy - D ^(mx ■ 4)^ (1) 

(xy) X 



The random- anisotropy model in the strong-anisotropy limit 



2 



where Sx is a three-component spin variable, unit vector describing the 

local (spatially uncorrelated) random anisotropy, and D the anisotropy strength. In 
amorphous alloys the distribution of Ux is isotropic, since, in the absence of crystalline 
order, there is no preferred direction. 

Random anisotropy is a relevant perturbation of the pure Heisenberg model, the 
crossover exponent being 0d = 0.412(3). Therefore, random-anisotropy systems 
show a behavior that is different from that observed in pure Heisenberg systems. Even 
though the critical behavior of the RAM has been investigated at length in the last 
thirty years (see [3] for a review) we do not have yet a satisfactory picture of its 
critical behaviour. The Imry-Ma argument forbids the presence of a low-temperature 
phase with nonvanishing magnetization for d < 4. However, this does not exclude the 
appearance of a glass transition with a low-temperature phase characterized by quasi- 
long-range order (QLRO), i.e., a phase in which correlation functions decay algebraically 
jnj. Functional renormalization-group calculations [HI E] indicate that QLRO may set 
in for small values of D, in agreement with a Landau-Ginzburg calculation [Hj of the 
equation of state for D —>■ 0. In the opposite limit D —>■ oo the RAM Hamiltonian 
reduces to that of an Ising spin glass with a correlated bond distribution. Indeed, for 
D ^ oo one can write = (^xUx, where cXa; = ±1 is an Ising spin, and obtain the 
Hamiltonian 



We call this model strong random-anisotropy model (SRAM) (We set J = 1 without 
loss of generality). Model (j2I) differs from the usual Ising spin glass model (ISGM) in 
the bond distribution: in Hamiltonian (j2)) the random variables jxy on different lattice 
links are correlated. For instance, one has rinii^y = 1/27, where the product is over 
the links belonging to a given plaquette and the average is taken with respect to the 
distribution of the vectors Ux- 

An interesting hypothesis, originally put forward in [Hj, is that the SRAM transition 
is in the same universality class as that of the ISGM. This conjecture looks very plausible, 
since the SRAM is nothing but an Ising model with local disorder and frustration. In 
some sense we can think of the SRAM and of the ISGM as two different versions of the 
same model: in the SRAM disorder is associated with lattice sites, while in the ISGM 
disorder is associated with lattice bonds. They are analogous to the site-diluted and 
bond-diluted Ising model (TUl E] , whose Hamiltonian is given by (j21) with jxy = VxVy 
(site dilution) and jxy = Txy (bond dilution), r being a random variable such that r = 1 
(r = 0) with probability p (resp. 1 — p). Note that the site-diluted model can also 
be interpreted as a bond-diluted model with a correlated bond distribution, exactly 
as is the case for the SRAM. Nonetheless, there is little doubt — though a satisfactory 
numerical check is still missing — that the two models belong to the same universality 
class. 

Note that the SRAM is less frustrated than the standard ISGM, since in the SRAM 
rinia:?/ — 1/27. This fact does not rule out the conjecture since it is known that maximal 




]xy — '^x ' 



(2) 
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frustration is not necessary to obtain glassy behaviour. For instance, the random-bond 
Ising model with j^y = +1 with probability p and jxy = —I with probability 1 — p has 
a glassy low-temperature phase for 0.233 ^ p ^ 0.767. 

The identity of the SRAM and ISGM universality classes was confirmed in two 
dimensions by a renormalization-group calculation using the large-cell method |13j . 
though in this case the critical point is at T = 0. In three dimensions instead, the 
critical behaviour of the SRAM has been controversial for a long time. While, for small 
D, numerical simulations [TH UHl CEl CZl CHI provided some evidence of the existence of 
a finite-temperature transition (though QLRO was never observed), in the SRAM even 
the existence of a finite-temperature transition was in doubt 1,8.^. 

In [in] we study again the SRAM and find good evidence for the existence of a finite- 
temperature transition. The corresponding critical behaviour turns out to be compatible 
with the conjecture of p. Note that at the transition only the overlap variables, which 
are the usual order parameters at a spin-glass transition, become critical. Magnetic 
quantities do not show any critical behaviour, though we expect nonanalyticities induced 
by the critical modes. Thus, on both sides of the transition the system is paramagnetic. 
It is unknown whether this paramagnetic phase survives up to T = 0. 

Note that our results predict an ISGM transition also in a generalized SRAM 
in which the vectors u are dimensional, for > 3. Indeed, as A^ increases, the 
bond correlation decreases and, for A^ — > oo, one reobtains the ISGM although with 
a nonstandard bond distribution. For N = 2 instead a ferromagnetic phase transition 
would be possible since the system is less frustrated than the case we have considered. 
Note that such a transition would only be observed in correlations of exC^, where the 
Ising variables ex should be chosen such as to have the couplings exeyjxy ferromagnetic on 
a maximal set of lattice links. This ferromagnetic transition would not violate the Imry- 
Ma argument, since order in the a^. variables does not imply order in the continuous 
variables Sx- Most likely, as for A^ = 3, the low-temperature phase would still be 
paramagnetic even if some Ising variables magnetize. 

A question that remains open is the behaviour of the RAM for finite anisotropy 
D. If there is indeed a low-temperature phase with QLRO for small D as predicted in 
jHllZllHj, then there should be a critical value D* such that ISGM behaviour is observed 
only for D > D*. Nothing is known about D* and we cannot even exclude that D* = oo, 
so that ISGM behaviour is observed only for model 

2. Results 

In Ref. ^2] we study the critical behaviour of the SRAM by means of Monte Carlo 
(MC) simulations. Since the model is essentially a spin glass we focus on the so-called 
overlap variables ^xTx, where ax and Tx are associated with two different replicas of the 
model with the same bond variables. For the SRAM one can also consider the standard 
magnetic variables Sx = CxUx- We find, in agreement with the results of that these 
quantities are not critical, i.e. on the low-temperature side of the transition the system 
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is still paramagnetic. We study the behaviour of the SRAM in the high-temperature 
phase. This reduces the algorithmic problems — the MC algorithm becomes very slow 
as temperature is reduced — and allows us to consider lattices of size up to L = 30. 
For the MC dynamics, we combine the Metropolis algorithm with the random-exchange 
(parallel-tempering) method [201 Ell- 
in order to verify whether the system has a critical behaviour, we looked for 
the occurence of finite-size scaling (FSS). For this purpose, we considered the ratios 
Ra = A{P, sL)/A{(3, L) with A = x, (both quantities are associated with the two- 
point function of the overlap variables), fixing s = 3/2. If FSS holds, as L increases 
all points should eventually fall on a universal curve depending on C,{j3,L)/L. Our 
results, reported in Figures [T[|21 show that this happens: The data corresponding to the 
pairs L = 16, 24 are only slightly different from those with L = 20, 30, the difference 
being much less than that observed by comparing data with L = 16, 24 and L = 8, 12. 
Therefore, the data strongly suggest that FSS holds, (though, for L < 30, scaling 
corrections are significant) and therefore, that the system becomes eventually critical. 
At the critical point R^{l3c,L) = s. Looking at Fig. |21we see that all MC data satisfy 
R^{I3,L) < 3/2 = s, which indicates that they all lie in the high-temperature phase. 
This allows us to set a lower bound on the position of the critical point, (3c ~ l-l- 

In order to determine the critical properties of the model, i.e. critical point and 
critical exponents, we used the iterative method which was introduced in |221 and 
generalized in [221 include scaling corrections. It allowed us to obtain infinite- volume 
estimates of x and ^ up to ^oo ~ 20 (^oo is the infinite- volume second-moment correlation 
length) in the high-temperature phase. The starting point is the FSS relation 

= Fa{s, m L)/L) + L-Ga{s, m L)/L), (3) 

valid for any long-distance quantity. Here s is an arbitrary number (in our calculations 
we fixed s = 3/2), Fa{s, z) and Gyi(s, z) are universal scaling functions, and is a to-be- 
determined correction-to-scaling exponent. By iterating (jHl) it is possible to extrapolate 
the finite- volume estimates A{(3,L), obtaining the infinite- volume value A^{(3). 

Our results show that ioo{.(3) increases quite rapidly in the range we have studied, 
0.80 ^ j3 ^ 1.0, confirming that the system eventually becomes critical. Fitting the 
infinite-volume results, we are able to determine the critical point and critical exponents 
associated with the critical behaviour of the overlap observables. We obtain: 

13c = 1.08(4) 
r]o = -0.24(4) 

lyo = 2.4(6) (4) 
7o = 5.3(1.3). 

The critical exponents are defined by x ~ and ^ ~ {f3c — ■ The suffix a is 

introduced to remind that they refer to the overlap variables and not to the magnetic 
ones. In the analysis it is crucial to include corrections to scaling in the FSS (corrections 
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Figure 1. The FSS curve of the susceptibility x Figure 2. The FSS curve of the correlation 
for s — 3/2. length ^ for s = 3/2. 
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Figure 3. The quartic cumulant Bg of the Figure 4. The ratio ^/L for several lattice sizes, 
overlap parameter for several lattice sizes. 



behave as and in fits of infinite- volume quantities (corrections behave as as 

^oo oo). The exponent u was determined by studying the critical behaviour of two 
universal ratios that involve the four-point and the two-point correlation function of the 
overlap variables. We obtained u = 1.0 ± 0.4. 

Estimates (jl} are close to those obtained for the ISGM and thus support the 
conjecture that the SRAM transition is in the same universality class as that of the 
ISGM. Ref. H quotes Uo = 2.22(15), rjo = -0.349(18), while |2S1 reports Uo = 2.39(5) 
and rjo = —0.395(17). Other estimates are reported in Table 1 of With the quoted 
error bars there is a small discrepancy between our estimate of t]o and those of fI^ ESI • 
This difference should not be taken too seriously, since there are similar discrepancies 
among the estimates obtained by different groups for the bimodal ISGM, see Table 1 of 

In MC simulations the critical point is often determined by considering the crossing 
point /Across of the Binder cumulant Bg = Jh/JIq^, where Hk = {(^2x1^)'')^ ratio 
^/L: indeed, /Across — >■ /3c as L — > oo. Our data for Bg and C,/L, which are reported in 
Figures El and m are compatible with /Across ~ 1.08. Indeed, the estimates of both Bg 
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and ^/L at fixed L get closer as (3 increases towards 1.08, although a crossing point is 
not clearly observed for those values of L (L < 16) that extend up to /3 = 1.10. This is 
probably due to scaling corrections that are particularly large in the SRAM. 

We can perform a more quantitative check by using the results of [22] for the 
critical-point values B* and (^/L)*. They quote: B* = 1.475(6) (bimodal distribution) 
and B* = 1.480(14) (Gaussian distribution); (^/L)* = 0.627(4) (bimodal distribution) 
and {^/L)* = 0.635(9) (Gaussian distribution). These results are compatible with ours 
for Bg and ^/L close to /? = 1.08. For /? = 1.07 we have Bg = 1.411(4) (L = 12), 
Bg = 1.434(6) (L = 16), and ^/L = 0.662(4) (L = 12), ^/L = 0.648(4) (L = 16). These 
results are very close to the ISGM estimates and show the correct trend. They are 
therefore consistent with the existence of a critical point at /5 = 1.08 ± 0.04 and with 
the conjecture that the ISGM and the SRAM belong to the same universality class. 
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